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Abstract 



We have performed ab initio mixed-states and sum-over-states calculations 
of parity nonconserving (PNC) electric dipole (El) transition amplitudes be- 
tween s-d electron states of Cs, Fr, Ba + , and Ra + . For the lower states of 
these atoms we have also calculated energies, El transition amplitudes, and 
lifetimes. We have shown that PNC El transition amplitudes between s-d 
states can be calculated to high accuracy. Contrary to the Cs 6s-7s transi- 
tion, in these transitions there are no strong cancelations between different 
terms in the sum-over-states approach. In fact, there is one dominating term 
which deviates from the sum by less than 20%. This term corresponds to an 
s-pi/2 weak matrix element, which can be calculated to better than 1%, and 
a Pi/2-^3/2 El transition amplitude, which can be measured. Also, the s-d 
amplitudes are about four times larger than the corresponding s-s transitions. 
We have shown that by using a hybrid mixed-states / sum-over-states approach 
the accuracy of the calculations of PNC s-d amplitudes could compete with 
that of Cs 6s-7s if Pi/2-^3/2 El amplitudes are measured to high accuracy. 



PACS: 32.80.Ys,31.15.Ar 
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I. INTRODUCTION 



Precise low-energy experiments on parity nonconservation (PNC) in atoms provide a test 
of the standard model of elementary particle interactions. By measuring PNC electric dipole 
(El) transition amplitudes, the value of the nuclear weak charge can be extracted by compar- 
ison with calculations. In a recent PNC experiment with cesium |IJ the PNC El transition 
amplitude between the 6s and 7s states has been determined with an unprecedented accu- 
racy of 0.35%. However, interpretation of the experiment is limited by the accuracy of the 
atomic calculations. Since 1989, calculations of the 6s-7s transition in Cs have been at the 
1% level At this level of accuracy the value of the nuclear weak charge is consistent with 
that predicted by the Standard Model. Recent measurements of values relevant to the PNC 
El amplitude (El transition amplitudes, hyperfme structure constants) are in much better 
agreement with the calculated values than they were ten years ago. From this it has been 
claimed that the accuracy of the calculated PNC El amplitude is 0.4% ||. Re-interpreting 
the Cs measurement with the higher accuracy, while using the calculations [0|§, the value 
of the nuclear weak charge gives a 2.5 o deviation from the Standard Model prediction [|J. 
However, inclusion of the Breit interaction into the calculations reduces the deviation by 
about 1 o H||. Note that these measurements give the best limits on new physics beyond 
the Standard Model, such as extra Z bosons, leptoquarks, composite fermions [0-0. 

One obviously needs an independent confirmation of the Cs result. In this paper we show 
that the accuracy of calculations of PNC El transitions between s-d states of Cs, Fr, Ba + , 
and Ra + could compete with that of the Cs 6s- 7s transition. The experiment for the 6s-5<i 
transition in Ba + is already in progress |]10|. 



II. MANY-BODY CALCULATIONS 

We perform calculations for iV-electron atoms with one external electron above closed 
shells. The calculations start from the relativistic Hartree-Fock (RHF) method in the V^ -1 
approximation. The single-electron RHF Hamiltonian is 

H = ca ■ p + (p - l)c 2 - Z/r + V N ~\ (1) 

a and j3 are Dirac matrices and p is the electron momentum. The accuracy of RHF energies 
is of the order of 10% for heavy atoms like Cs, Fr, Ba + , and Ra + . 

In order to obtain more realistic wavef unctions, electron-electron correlations must be 
taken into account. Correlation corrections to the electron orbitals are calculated using 
the "correlation potential" method |TT[. This method corresponds to adding a non-local 
correlation potential £ to the potential V^" 1 in the RHF equation (|]) and then solving for 
the states of the external electron. The correlation potential is defined such that its average 
value coincides with the correlation correction to energy, SEi = (^jlE^). The correlation 
potential is calculated by means of many-body perturbation theory in the residual Coulomb 
interaction 

N N -i N 

U = H-Y. = E T-^-l ~ E V*- 1 ^), (2) 

i=l i<j I 1 * l j\ i=l 
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where H is the exact Hamiltonian of an atom. The lowest-order correlation diagrams 
(second-order in U) are presented in Fig. |I|. At this level of calculation the accuracy for 
energy levels is about 1%. 

Using the correlation potential method and the Feynman diagram technique we include 
three series of higher order diagrams which are calculated in all orders of perturbation 
theory [^,|TJ,|T^]. These are screening of the electron-electron interaction, the hole-particle 
interaction, and chaining of the self-energy operator E. The electron-electron screening 
(see Fig. |j) and the hole-particle interaction (Fig. ||]) are incorporated into the self-energy 
operator E (Fig. [|). Chaining of the self-energy operator to all orders (Fig. |5]) is then 
calculated by adding E to the Hartree-Fock potential V N ~ l and solving the equation 

(H + E - e)V = (3) 

iteratively for the states of the external electron. In this way "Brueckner" energies and 
orbitals are obtained. These energies have an accuracy of the order of 0.1%. The wavefunc- 
tions can be further modified by placing a coefficient before E such that the corresponding 
energy coincides with the experimental value. This fitting of the Brueckner orbitals can be 
considered as a way of including higher-order diagrams into the calculations. 

We use the time-dependent Hartree-Fock method (which is equivalent to the random- 
phase approximation with exchange) to calculate the interaction of external fields with 
atomic electrons. In this paper we deal with two external fields: the electric field of the 
photon (El transition amplitudes) and the weak field of the nucleus. In the RHF approxima- 
tion the interaction between an external field H ext and atomic electrons is (ipi F \H ext \ip2 F )^ 
where ipi F and ipi are RHF orbitals. Inclusion of the polarization of the atomic core 
by an external field is reduced to the addition of a correction 5V (which is the correction 
to the Hartree-Fock potential due to the interaction between the core and the external 
field) to the operator which describes the interaction, (ip^ F \H cxt + 5V\ip2 F ). To include 
"Brueckner-type" correlation corrections the RHF orbitals are simply replaced by Brueckner 
ones, (ipi r | H ext + 8V\ip2 V )- The Brueckner-type correlations give the dominant corrections 
to the RHF approximation. They correspond to diagrams in which the interactions occur 
in the external lines of the self-energy operator (see, e.g., Fig. ||). Those diagrams in which 
the El interaction occurs in the internal lines are known as "structural radiation", while 
those in which the weak interaction occurs in the internal lines are known as the "weak 
correlation potential" (see, e.g., Fig. |^). There is also a correction to the amplitudes arising 
from the normalization of states [ITf. The structural radiation, weak correlation potential, 



and normalization contributions are suppressed by the small parameter E ext /E cove ~ 1/10, 
where E ext and E- int are excitation energies of the external and core electrons, respectively. 
The nuclear spin-independent weak interaction of an electron with the nucleus is 

G 

Hw = -^^p(r)Q w -f 5 (4) 

where Gf is the Fermi constant, Qw is the weak charge of the nucleus, 75 is a Dirac matrix, 
and p(r) is the density of the nucleus. Parity nonconserving El transition amplitudes, arising 
due to the simultaneous interaction of atomic electrons with the nuclear weak charge and 
the photon field, can be calculated using two methods: from a mixed-states approach; or 
from a sum-over-states approach, in which experimental values (energies and El transition 
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amplitudes) can be explicitly included. Contributions to PNC El transition amplitudes are 
presented diagrammatically in Figs. §, 0. 

In the mixed-states approach the PNC El transition amplitude between the states ns 
and (n — l)d 3 / 2 , n = 6 for Cs and Ba + , n = 7 for Fr and Ra + , is given by 

Elp N c = (i>(n-X)d\HEl + SV E l\Sll) ns ) + (lf)( n -.i) d \H W + SVw\i>ns) + (i>(n-l)d\8V E lw\i>ns) , (5) 

where 5ip and 5Vw are corrections to single-electron wavefunctions and the Hartree-Fock 
potential caused by the weak interaction, ip and 5V E \ are corrections to wavefunctions and 
the Hartree-Fock potential caused by the electric field of the photon, 5V E rw is the correction 
to the core potential due to the simultaneous action of the weak field and the electric field 
of the photon; the wavefunctions tp( n -i)d an d ipns correspond to Brueckner orbitals, and the 
corrections 5ip and ip are found by solving the equations: 

(H + t - e)5i[> = -{H w + 8V W )^ 
{H + £-e)<iP = -{H E1 + 6V E1 )ip. 

This method is equivalent to calculating the diagrams presented in Fig. ^| (with £ chained 
to all orders, Fig. |5|) with the inclusion of the core polarization diagrams presented in Fig. || 
Parity nonconserving El transition amplitudes between the states ns and (n— l)d in the 
sum-over-states approach have the form 

ri v-^ ((n - l)d 3 / 2 \H E i + SVeiWp!^) (n'p 1/2 \H w + 5V w \ns) 
uipNc - 2^ p _ p 

„ ((n - l)d 3/2 \H w + SV w \n'p 3 / 2 ) (n'p 3/2 \H E i + 5V E1 \ns) 

2-1 Tp _ rp > \°) 

where the sum is taken over a complete set of pi/ 2 and p 3 / 2 states. 

Note that the sum-over-states approach should also include the states with double- 
excitations like, for example in Cs, (5p 6 ds\Hw\5p &p7 s) . In the mixed-states approach these 
states are included, for example, in the last term of Eq. || (see also diagram (c) of Fig. ||). 
These exotic states contribute due to their mixing with the single-excited electron states. 
This means that the mixed-states calculation @ is more complete than the sum-over-states 
([]) unless the high-energy states with two or more excited electrons are included into the 
sum. 

However, the accuracy of pure ab initio calculations for s-d transitions is not very good 
because of the huge correlations for d-states. On the other hand, we will see in the next 
section that this problem can be avoided in the sum-over-states approach by using exper- 
imental values for the p-d El transition amplitudes. Therefore, the best accuracy can be 
achieved when both methods are combined. Substitution of experimental values into the 
sum-over-states approach leads to a correction to the PNC amplitude which can be added 
to the mixed-states result. Following this procedure, it is possible to determine the PNC s-d 
amplitudes with an accuracy of about 1% (see discussion in the end of the next section). 
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III. RESULTS 



Hartree-Fock energies for Cs, Fr, Ba + , and Ra + are presented in Table |. These have an 
accuracy of the order of 10%. The Brueckner energies, including the three series of higher 
order diagrams, are also presented in Table |I[ These energies have an accuracy of the order 
of 0.1%. 

Electric dipole transition amplitudes between the states m and m' are calculated in length 
form, (m'\\H E1 \\m) = (m'\\r\\m) = C mm 'R, where C mm i are angular coefficients and R is 
the radial integral. In Table [II] we present radial integrals relevant to the sum-over-states 
calculation for Cs, Fr, Ba + , and Ra + . In this table we present the values obtained in the 
RHF approximation and show the contribution of core polarization to the RHF integrals; we 
also present the (unfitted) Brueckner results and the contributions arising from structural 
radiation and normalization of states. In Table |TJ we present radial integrals between 
the lower states of the four atoms which are calculated with fitted Brueckner orbitals and 
with structural radiation and normalization contributions included. Experimental values 
for Cs and Fr are presented in Table |V|. The Cs transitions 5rf 3/ / 2 -6p3/2 and 5<i3/2-6pi/2 
were extracted from the measurement [56[] of the 5g?3/2 lifetime, r = 909(15) ns by assuming 



that the ratio of the calculated radial integrals corresponds to the ratio of the experimental 
values. This assumption was also used to obtain the Fr 7d 3 / 2 -7p 3 /2 and 7d 3 / 2 -7pi/2 radial 



integrals from the measured lifetime r = 73.6(3) ns [28|. With the exception of the Cs 6s 



7pi/2 transition, the calculations of s-pi/2 radial integrals agree with experiment at the level 
of 0.1%. The poor accuracy of the 6s-7pi/2 radial integral is due to the fact that the main 
RHF contribution is very small and the relative contribution of all corrections is large. The 
5d-6p radial integrals for Cs have poor accuracy, deviating from experiment by about 4%. 
This is indicative of the poor calculation of d-states due to very large correlation corrections. 
The accuracy for the Fr 7d-7p radial integrals is about 1%. The reason that this accuracy is 
better than that for Cs 5d-6p is because the accuracy of the higher <i-levels (here 7d rather 
than 6c?) is better due to smaller correlation corrections. 

We have calculated the lifetimes of the low-lying states of Ba + and Ra + . The nd^/2 states 
of Ba + (n = 6) and Ra + (n = 7) decay directly to the ground state via the E2 transition; the 
nd^/2 states decay via both the E2 and Ml transitions. Lifetimes of the Ba + 5g? states are 
presented in Table M. The calculations were performed with fitted Brueckner orbitals; core 
polarization, structural radiation, and normalization contributions were included into the E2 
transition amplitudes. The calculations are in good agreement with experiment. We have 



also presented the calculations performed by Guet and Johnson [29]. It appears that for the 



state 5g?5/2 they have not taken into account the Ml transition. This seems to be the reason 
for the discrepancy between the lifetime calculations for this state. From our calculations 
it is seen that inclusion of the Ml transition effectively decreases the lifetime of the bd^/ 2 
state from 36.3 s to 30.3 s. For Ra + we obtained the lifetimes r = 0.641 s and r = 0.302 s 
for the states 6g? 3 / 2 and 6g? 5 / 2 , respectively; these were calculated in the same way as for the 
Ba + 5d states. The lifetimes for states of Ra + have not been measured. Lifetimes of all 
other states are strongly dominated by El transitions and so can be calculated using the 
radial integrals presented in Table |TJ . The calculated lifetimes of the Qp states of Ba + are in 
excellent agreement with experiment and calculations by Guet and Johnson (see Table [V]). 
The mixed-states results for the El PNC transition amplitudes are listed in Table [VT|. 
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The results of the sum-over-states calculation, and the contributions of the six terms corre- 
sponding to the summation of the n-(n + 2) p states, are presented in Table [VII . In both 
calculations the contributions of structural radiation, weak correlation potential, and nor- 
malization of states are not included. The s-d PNC amplitudes are up to about four times 
as large as their corresponding s-s amplitudes. Furthermore, unlike the contributions to the 
sum-over-states calculation in Cs 6s-7s, in which the dominant contribution is about twice 
as large as the final result due to strong cancelations between three major terms in the sum, 
the PNC s-d transitions in Cs, Fr, Ba + , and Ra + are strongly dominated by a single term. 
In each case this term corresponds to ((n — 1)d^/2\H E i\npi/ 2 ) (npi/ 2 \Hw\ ns ) I '(E ns — E npi/2 ); 
this term is different from the sum by less than 20%. 

Because the Cs Qpi/2-5d3/2 El transition and energies are known we can correct the 
mixed-states PNC result. Replacing the calculated values by these experimental values in 
the dominating term of the sum (^j) for Cs increases this term (and the total sum) by about 
4%. This correction is mostly due to the difference between the calculated and experimental 
El amplitude (see Table [TV]). From the 1% accuracy of calculations of hyperfine structure 
constants for s and pi/2 states JIB we can expect that the accuracy of the s-pi/2 weak matrix 
elements in this calculation is also about 1%. Therefore, we can say that the uncertainty in 
the calculated ab initio s-d E1 PNC amplitudes is dominated by the uncertainty of the p-d 
El matrix elements and constitutes about 4% for Cs and about the same value, or a little 
more, for other atoms. 

The correction to the 6s-5d Elp^c amplitude in Cs discussed in the previous paragraph 
is 0.126 ieciB{—Qw /N). When it is added to the mixed-states result, the new value is 
ElpNc(6s — 5d) = 3.75 ieap{—Qw /N). Since using the experimental 6p-5d El amplitude 
removes the main source of uncertainty, the accuracy of the modified result must be con- 
siderably better than 4%. Assuming high accuracy of p-d transition amplitudes, one can 
say that the uncertainty is now dominated by the uncertainty of calculated s-p weak matrix 
elements which is about 1% (however, we believe that this accuracy can be improved beyond 
1% with the inclusion of weak correlation potential and normalization contributions). More 
rigorous calculations and a more detailed analysis of the accuracy will be carried out when 
the need arises from the progress in experiments. We expect that calculations of Elp^c 
amplitudes for Cs, Fr, Ba + , and Ra + , with empirical corrections, can reach an accuracy of 
about 1%. 



IV. CONCLUSION 

We have calculated the PNC El transition amplitudes between s-d states of Cs, Fr, Ba + , 
and Ra + . Generally, high accuracy cannot be reached in purely ab initio calculations of these 
transitions due to the poor accuracy of d-states. However, we have shown from a sum-over- 
states calculation that, unlike the Cs Qs-7s transition, the s-d transitions we have mentioned 
are strongly dominated by a single term in the sum. Moreover, this term corresponds to an 
s-pi/2 weak matrix element, which can be calculated with an accuracy of better than 1%, and 
a ^1/2-^3/2 El transition amplitude, which can be taken from experiment. The need to reach 
high accuracy for d states is therefore avoided. In addition to this, PNC s-d transitions are 
larger than the corresponding s-s transitions. The mixed-states calculation can be modified 
by correcting the terms in the sum-over-states by inserting experimental El transitions and 
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energies. If pi/2-^3/2 El transition amplitudes are measured to high accuracy, we believe 
that the accuracy of the calculations of s-d PNC transitions for Cs, Fr, Ba + , and Ra + can 
reach 1%. 
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TABLES 



TABLE I. Energy levels (ionization potentials) of the lower states of Cs, Fr, Ba + and Ra + in 



units —cm 









Cs 






Ba+ 




State 


RHF 


Brueckner 


Experiment a 


RHF 


Brueckner 


Experiment a 


6s 


27954 


31420 


31407 


75339 


80813 


80687 


7s 


12112 


12851 


12871 


36852 


38333 


38332 


8s 


6793 


7082 


7090 


22023 


22651 


22662 


6pi/ 2 


18 


790 


20275 


20228 


57265 


60581 


60425 


7pi/ 2 


9223 


9643 


9641 


30240 


31332 


31297 


8P1/2 


5513 


5701 


5698 


18848 


19378 


19351 


6P3/2 


18 


389 


19708 


19674 


55873 


58860 


58734 


7P3/2 


9079 


9460 


9460 


29699 


30704 


30676 


8P3/2 


5446 


5618 


5615 


18580 


19075 


19051 


54/2 


14138 


17023 


16907 


68139 


76402 


75813 


64/2 


7920 


8824 


8818 


33266 


34740 


34737 


74/2 


4965 


5362 


5359 


20251 


20871 


20887 


54/2 


14163 


16915 


16810 


67665 


75525 


75012 


64/2 


7921 


8781 


8775 


33093 


34536 


34532 


74/2 


4963 


5341 


5338 


20167 


20777 


20792 








Fr 






Ra+ 




State 


RHF 


Brueckner 


Experiment b 


RHF 


Brueckner 


Experiment a 


7s 


28768 


32841 


32849 


75900 


81960 


81842 


8s 


12282 


13071 


13116 


36861 


38405 


38437 


9s 


6858 


7164 


7178 


22004 


22659 


22677 


7P1/2 


18856 


20674 


20612 


56879 


60681 


60491 


8P1/2 


9240 


9730 


9736 


30053 


31244 


31236 


9pi/2 


5521 


5737 




18748 


19332 




7^3/2 


17656 






ozyuo 


OO 1 04 


OOUOO 


8^3/2 


8811 


9180 


9191 


28502 


29447 


29450 


9ft/2 


5319 


5486 




17975 


18462 


18432 


64/2 


13826 


16610 




62356 


70149 


69758 


74/2 


7725 


8583 


8604 


31575 


33060 


33098 


84/2 


4857 


5241 


5248 


19451 


20079 


20107 


64/2 


13925 


16413 




61592 


68449 


68099 


74/2 


7747 


8496 


8516 


31204 


32569 


32602 


84/2 


4863 


5197 


5203 


19261 


19849 


19868 


a Taken from 


14 


I- 



'Measured in Refs. |§-§1|. 
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TABLE II. Calculated radial integrals (a.u.) for Cs, Fr, Ba + , Ra + . We present RHF values, 
RHF with core polarization, the Brueckner result with core polarization included, and structural 
radiation and normalization of states; 0.0 signifies that the value is smaller than the number of 
figures specified. 



Atom 


Transition 


RHF 


RHF + 


Brueckner + 


Structural 


Normal- 








core polar- 


core polar- 


radiation 


ization 








ization 


ization 




of states 


Cs 


6s 1/2 - 6p 3/2 


-6.432 


-6.074 


-5.500 


-0.028 


0.047 




6s 1/2 - 7f> 3 /2 


-0.602 


-0.440 


-0.463 


-0.013 


0.003 




6s 1/2 - 8f> 3/2 


-0.245 


-0.143 


-0.162 


-0.008 


0.001 




6pi /2 - 5d 3 / 2 


7.775 


7.481 


6.050 


0.026 


-0.050 




iTli In F»/7o In 

'Fl/2 ou 3/2 


-3.498 


-3.591 


-1.742 


0.011 


0.011 




8P1/2 - 5d 3 /2 


-0.860 


-0.916 


-0.556 


0.007 


0.003 


Fr 


7s 1/2 - 7p 3 /2 


-6.140 


-5.739 


-5.128 


-0.032 


0.051 




7s 1/2 - 8p 3/2 


-0.949 


-0.760 


-0.748 


-0.015 


0.006 




7s 1/2 - 9?3/2 


-0.452 


-0.332 


-0.336 


-0.009 


0.003 




7p 1/2 - 6(f 3 / 2 


7.986 


7.613 


6.256 


0.035 


-0.061 




8?}-i In 6f/o In 


-4.005 


-4.116 


-2.249 


0.014 


0.015 




9pi/ 2 - 6d 3/2 


-0.941 


-1.008 


-0.704 


0.008 


0.004 


Ba+ 


6si /2 - 6p 3/2 


-4.744 


-4.314 


-4.056 


-0.032 


0.046 




6si/ 2 - 7j> 3 /2 


-0.226 


-0.036 


-0.015 


-0.013 


0.0 




6si/ 2 - 8p 3/2 


-0.068 


0.054 


0.076 


-0.007 


-0.001 




6pi/ 2 - 5d 3 / 2 


3.244 


2.964 


2.634 


0.026 


-0.034 




7pi /2 - 5d 3/2 


0.304 


0.184 


0.225 


0.010 


-0.002 




8pi /2 - 5d 3/2 


0.169 


0.092 


0.099 


0.007 


-0.001 


Ra+ 


7s 1/2 - 7p 3/2 


-4.624 


-4.154 


-3.885 


-0.035 


0.048 




7s 1/2 - 8p 3/2 


-0.541 


-0.320 


-0.286 


-0.015 


0.003 




7s 1/2 - 9p 3/2 


-0.243 


-0.098 


-0.067 


-0.009 


0.001 




7pi/ 2 - 6d 3/2 


3.851 


3.448 


3.067 


0.037 


-0.047 




8pi /2 - 6d 3/2 


0.091 


-0.075 


0.009 


0.014 


0.0 




9pi /2 - 6d 3 / 2 


0.075 


-0.030 


-0.014 


0.008 


0.0 
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TABLE III. Radial integrals (a.u.) for states of Cs, Fr, Ba + and Ra + . Fitted Brueckner 



orbitals are used; core polarization, and structural radiation and normalization of states are also 
included. 







6P1/2 


7^1/2 


8pi/ 2 


6^3/2 


7p 3 /2 


8^3/2 


Cs 


6s 


-5.508 


-0.313 


-0.081 


-5.482 


-0.471 


-0.171 




7s 


5.211 


-12.605 


-1.137 


5.625 


-12.383 


-1.419 




8s 


1.266 


11.386 


-21.753 


1.273 


12.163 


-21.252 




54/2 


6.072 


-1.785 


-0.560 


6.120 


-1.579 


-0.506 




64/2 


-3.696 


15.558 


-4.234 


-4.082 


15.612 


-3.726 




74/2 


-1.806 


-5.632 


27.670 


-1.918 


-6.379 


27.773 




54/2 


- 


- 


- 


6.190 


-1.641 


-0.522 




64/2 


- 


- 


- 


-3.986 


15.699 


-3.886 




74/2 


- 


- 


- 


-1.893 


-6.186 


27.876 


Ba+ 


6s 


-4.054 


0.121 


0.141 


-4.048 


-0.030 


0.063 




rj 

IS 


3.053 


o COO 

-8.583 


-0.139 


3.362 


o AHA 

-8.464 


-0.401 




OS 


U.oDo 


« nsn 

O.UoU 


i /i i ^ 

-14.10O 


U.ooo 


0.004 


1 1 SS/I 
-lo.oo4 




^A 

54/2 


O RAR 


n oof; 

U.ZZO 


n 1 HQ 


Z.0o4 


n ooc; 
U.zoO 


U.loO 




64/2 


A 09/1 

-4.zo4 


7 A CO 

7 .488 


0.101 


A con 
-4.520 


1 .311 


U.ZOZ 




74/2 


-1.189 


- / .134 


1 O AC\A 

lo.4y4 


1 1 7n 
-1.17U 


7 R1 A 

-l.biA 


-i o i nc 
13.195 




KA 

54/2 








O RKQ 


n 070 
yj.ziv 


U.loo 




6d5/2 








A ARQ 
-4.40t* 


7 A 1 S 
( .410 


n oik 
U.zoO 




7A 

'4/2 








-1.1800 


- ( .oo ( u 


lo.ooZ 






7P1/2 


8^1/2 


9P1/2 


7p 3 /2 


8^3/2 


9^3/2 


Fr 


7s 


-5.242 


-0.331 


-0.093 


-5.107 


-0.748 


-0.343 




8s 


5.217 


-12.326 


-1.206 


6.484 


-11.536 


-1.947 




9s 


1.256 


11.428 


-21.382 


1.215 


13.777 


-19.660 




64/2 


6.237 


-2.229 


-0.691 


6.417 


-1.553 


-0.516 




74/2 


-3.014 


15.893 


-5.488 


-4.186 


16.175 


-3.829 




84/2 


-1.606 


-4.296 


27.952 


-1.975 


-6.546 


28.464 




64/2 


- 


- 


- 


6.576 


-1.684 


-0.549 




74/2 






_ 


-3.974 


16.368 


-4.184 




8 4/2 


- 


- 


- 


-1.917 


-6.106 


28.695 


Ra+ 


7s 


-3.948 


0.108 


0.142 


-3.877 


-0.294 


-0.082 




8s 


3.104 


-8.523 


-0.168 


4.038 


-8.071 


-0.891 




9s 


0.867 


6.167 


-14.122 


0.897 


7.825 


-13.137 




64/2 


3.074 


0.011 


-0.011 


2.913 


0.245 


0.110 




74/2 


-3.774 


8.262 


-0.412 


-4.662 


7.801 


0.256 




84/2 


-1.240 


-6.152 


14.649 


-1.242 


-7.812 


13.854 




64/2 








3.109 


0.224 


0.105 




74/2 








-4.515 


8.087 


0.117 




84/2 








-1.281 


-7.476 


14.276 
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TABLE IV. Calculated (Table |T|) and experimental radial integrals (a.u.). 



Atom 


Transition 


Calc. 


Exp. 


Cs 


6s-6jDi in 


-5.508 


-5.497(8) a 




6s-6jDq/o 

Jro/ Z 


-5.482 


-5.476(6) a 




6s-7pi /o 


-0.313 


-0.348(3) b 




7«S-6?9i /o 


5.211 


5.185(27) c 




7s-6pq/o 

i O / Z 


5.625 


5.611(27) c 




7s-7pi /o 

Ir 1 / Z 


-12.605 


-12.625(18) d 






-12.383 


-12.401(17) d 




5(iq/9-6pi /o 

O / Z .f -L / Z 


6.072 


6.31(5) e 




5d 3 / 2 -6p 3 /2 


6.120 


6.36(5) c 




54/2-6^3/2 


6.190 


6.40(2) e 


Fr 


7s-7p 1/2 


-5.242 


-5.238(10) f 




7s-7p 3/2 


-5.107 


-5.108(13) f 




7d 3 / 2 -7pi/ 2 


-3.014 


-3.05(1) « 




7ci 3 /2-7p 3 /2 


-4.186 


-4.24(2) « 




7d 5 / 2 -7p 3 /2 


-3.974 


-4.02(8) 8 



a Reference 22|| . 
b Reference [23|. 
c Reference p4 |. 
d Reference [25|. 
e Reference |p6 |. 
f Reference [27]. 
g Reference 281. 



TABLE V. Lifetimes of low-lying states of Ba^ 



Atom 



State 



T 



this work 



r cxpt 



Ba 4 



6P1/2 
6P3/2 



5d 
5d 



3/2 
5/2 



7.89 ns 
6.30 ns 
81.5 s 
30.3 s 



7.99 ns 
6.39 ns 
83.7 s 
37.2 s 



7.90(10) ns b 
6.32(10) ns b 
79.8(4.6) s c 
34.5(3.5) s d 



a Reference |2S] 
b Reference [30 
c Reference pl | 
d Reference [32 
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TABLE VI. Mixed-states results for PNC El transition amplitudes between the states 
ns-{n — l)d, ((n — l)d\El z \ns) , where n = 6 for Cs and Ba + and n = 7 for Fr and Ra + ; units are 
10- u iea B (-Q w /N). 





Cs 


Fr 


Ba+ 


Ra+ 


Mixed-states Elp^c 


3.62 


57.1 


2.17 


42.9 



TABLE VII. Results of the sum-over-states calculations for the PNC El transition amplitudes 
for Cs, Fr, Ba + , and Ra + . We present the contributions of the terms in the sum corresponding 
to the intermediate n-(n + 2)p states, the contribution due to all other intermediate p-states, and 
the total value; units 10~ n ieaB(— Qw/N)] 0.0 means that the term is smaller than the number of 
figures specified. 







(d\ H E i 1 npi/2 > ("Pi/2 1 H w I s) 
E s -E nPl/2 




{d\Hw\np 3 /2)(np- s/2 \HEi\s) 
E d-E np . j/2 


Cs 


n=6 


3.154 




0.728 




n=7 


-0.258 




-0.013 




n=8 


-0.047 




-0.002 




Other 




0.197 






Total 




3.76 




Fr 


n=7 


59.78 




5.19 




n=8 


-6.13 




-0.15 




n=9 


-1.10 




-0.03 




Other 




1.95 






Total 




59.5 




Ba+ 


n=6 


2.036 




-0.264 




n=7 


0.045 




-0.001 




n=8 


0.012 




0.0 




Other 




0.511 






Total 




2.34 




Ra+ 


n=7 


40.69 




-2.33 




n=8 


0.11 




-0.05 




n=9 


0.02 




-0.01 




Other 




7.47 






Total 




45.9 
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FIGURES 



FIG. 1. Second-order correlation diagrams for the valence electron (S operator). Dashed line 
is the Coulomb interaction. Loop is the polarization of the atomic core. 
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FIG. 2. Screening of the Coulomb interaction. 




FIG. 3. Hole-particle interaction in the polarization operator. 



FIG. 4. The electron self-energy operator with screening and hole-particle interaction included. 
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FIG. 5. Chaining of the self-energy operator. 
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FIG. 6. Brueckner-type correlation corrections to the PNC El transition amplitude; the crosses 
denote the weak interaction and the dashed lines denote the electromagnetic interaction. 
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(a) 
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(c) 5Z 

FIG. 7. Small corrections to the PNC El transition amplitude: external field inside the correla- 
tion potential. In diagrams (a) the weak interaction is inside the correlation potential (5T, denotes 
the change in S due to the weak interaction); this is known as the weak correlation potential. 
Diagrams (b,c) represent structural radiation (photon field inside the correlation potential). In 
diagram (b) the weak interaction occurs in the external lines; in diagram (c) both the weak and 
electromagnetic interactions occur in the internal lines. 



— * 1 7 X ; ; 1 ; 7-^ 1 

/ / / 

/ / / 

(a) (b) (c) 

FIG. 8. Examples of diagrams representing the polarization of the atomic core by external 
fields. (The diagrams we have presented are exchange diagrams; there are also direct diagrams.) In 
diagrams (a) and (b) the core is polarized by a single field (the dashed line denotes the El interaction 
and the cross denotes the weak interaction). Diagram (c) corresponds to the polarization of the 
core by both fields. 
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